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Abstract 

A matrix model is constructed which describes a chiral version of the large U (N) 
gauge theory on a two-dimensional sphere of area A. This theory has three separate 
phases. The large area phase describes the associated chiral string theory. An exact 
expression for the free energy in the large area phase is used to derive a remarkably 
simple formula for the number of topologically inequivalent covering maps of a sphere 
with fixed branch points and degree n. 
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1 Introduction 

Although QCD is a relatively old subject, the nonperturbative calculation of physically sig- 
nificant quantities in this theory is poorly understood. Two-dimensional pure QCD, however, 
does admit a nonperturbative exact solution [1-3] . Recent work [4-8] has shown further that 
there is string theory description of pure two-dimensional gauge theories with various gauge 
groups. It is hoped that some modification of this string picture will provide a systematic 
nonperturbative description of QCD in any dimension. 

A particularly interesting feature of two-dimensional QCD is that when the theory is 
considered on a two-sphere, there is a phase transition between the large and small area 
limits of the large gauge theory. It was shown by Douglas and Kazakov [^ that the string 
expansion is asymptotic to the full nonperturbative QCD2 partition function only when the 
dimensionless area is greater than vr^; below this value, the theory simplifies dramatically 
and is essentially a free theory described by a master field satisfying the Wigner semicircle 
law. 

The existence of this phase transition raises the question of whether an analogous tran- 
sition will occur in the four dimensional theory. Such a phase transition would mean that a 
calculation in the strong coupling regime would be unable to make useful predictions about 
physics at weak coupling, so it is certainly important to understand the nature of this type 
of phase transition. 



Gross and Matytsin [^ have studied the DK phase transition from the point of view 
of the weak coupling expansion. They found that for weak coupling (in two dimensions), 
there are nonperturbative instanton contributions to the partition function which vanish in 
the asymptotic expansion below the critical point but which give rise to an infinite instan- 
ton density above that point, driving the DK phase transition. An analogous calculation 
indicated that this mechanism might not cause a phase transition in four dimensions. 

The phase transition has also been studied from the string point of view in |Tl| . Analytic 



and numerical results indicate that the string expansion has a radius of convergence which 
coincides with the critical point. In the string language, the transition is driven by the 
entropy of branch point singularities in the string maps. (The formula for this entropy 



which was conjectured in [^ is proven in this paper.) 

Although these studies have given us some understanding of the reasons for the phase 
transition in two-dimensional QCD, there are still many unanswered questions. Studying 
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similar phase transitions in closely related theories may give insights into the structure of 
these transitions. 

An interesting feature of the large QCD theory in two dimensions is that the theory 
almost exactly factorizes into two copies of a simpler, "chiral" theory @, 0]. The chiral theory 
encapsulates most of the important features in the geometry of the underlying string theory, 
and seems worthy of study as a theory in its own right. 

The goal of this paper is to study a chiral sector of the two-dimensional gauge theory 
on the sphere and to determine its large N phase structure. We find that the chiral theory 
contains 3 distinct phases, as opposed to the 2 phases of the combined theory. Furthermore, 
we find that the same analytic formulae describe both the large and small area limits; in 



analogy with what was found in [|1T|. These formulae also reproduce the asymptotic string 
expansion. This suggests that the chiral theory, which is not yet well understood from the 
point of view of the large gauge theory, may have unexpected duality properties. 

In Section 2 we briefly describe the QCD2 string and review the large A^ theory and the 
Douglas-Kazakov phase transition. We define the chiral theory and give an Ansatz for the 
matrix model corresponding to this theory. Section 3 contains a saddle point evaluation of 
the chiral partition function. We find a class of single-cut saddle point solutions which agree 
order by order with the string expansion and which are defined for large and small areas, 
but not in an intermediate region. In Section 4, the single-cut solution is used to derive 
a previously conjectured formula for the number of sphere to sphere maps with n branch 
points. In Section 5 we find double-cut solutions for the saddle point equations and show 
that these solutions interpolate exactly between the regions of weak and strong coupling 
given by the single-cut solutions. Section 6 contains a discussion of our results and related 
issues. 



2 Preliminaries 

The partition function of two-dimensional Yang-Mills theory on a compact manifold with 
genus G and area A is given by the exact formula p|, |[ 

ZiG,A,N) = j[VA^]e-^^^^^^''*^^ 

= ^(dimi?)2-2^e-wC2(«), (2.1) 

R 
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where the sum is taken over all irreducible representations of the gauge group, with dimi? 
and C2{R) being the dimension and quadratic Casimir of the representation R. We absorb 
the coupling constant A = g'^N into the dimensionless area A throughout this paper. 

We are interested in studying the Yang-Mills theories with gauge groups U{N) and 
SU{N). The irreducible representations of these gauge groups are associated with Young 
diagrams with row lengths rij, where for U{N) the row lengths satisfy 

oo > ni > ■ • • > njy > — oo (2.2) 

and for SU{N) 

oo > rii > ■ ■ ■ > riN = 0. (2.3) 

We denote the sum of the row lengths hj n = J^f^i- In terms of the row lengths, the 
dimension of a representation R is given by 

dimi? = nfl-^^^^V (2.4) 



i>j 



For U{N) the quadratic Casimir is given by 

C2{R) = Nn + Y.ni{n,-2t+l), (2.5) 

i 

and for SU{N) we have 

C2{R) = Nn + Y.ni{m-2i + l)-—. (2.6) 

i 

In the asymptotic expansion in of the SU{N) partition function ( p.l|) was 
described in terms of a string theory of covering maps. This asymptotic expansion was 
constructed by summing over a subset of representations called "composite" large N repre- 
sentations; these are essentially the representations whose quadratic Casimir has a leading 
term of order A^. Composite representations are formed by taking the tensor product of a 
representation R corresponding to a Young diagram with a finite number of boxes, and a 
representation S corresponding to the complex conjugate of a representation S associated 
with another diagram with a finite number of boxes. The resulting large N expansion of 
( p.l|) essentially factorizes into two copies of a simpler, "chiral" string expansion. The chiral 
expansion arises from restricting the sum to only representations corresponding to diagrams 
with a finite number of boxes. The string geometry of the chiral theory is simpler than that 
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of the full theory because in the chiral theory all maps from the string world sheet to the 
target space have the same relative orientation, so the theory can be viewed as a theory of 
orientation-preserving string maps. The chiral QCD2 string theory has recently been con- 
structed by Cordes, Moore, and Ramgoolam in terms of a topological string theory [|12 



related work appears in [|T3|. 

In the string expansion of ( |2.1| ), a term of order A^^~^^ in the partition function is inter- 
preted as arising from a (possibly disconnected) string world sheet of genus g. The dimension 
of a representation associated with a Young diagram containing a finite number n of boxes 
has a leading term of order iV". Thus, when the genus G is greater than 1, the leading 
order nontrivial term in the chiral partition function is of order N"^'"^^ and arises from the 
fundamental representation (with Young diagram □). In the full SU{N) theory, both the 
fundamental representation and its conjugate contribute terms of order A^^~^*^. Because in 
these cases only a finite number of terms contribute to the leading order term in the partition 
function, the large N theory with G > 1 clearly cannot have a phase transition at leading 
order. When G = 1, the contribution of the dimension factor to the partition function 
vanishes and so all representations contribute to the A^'' term. However, in this case the 
contribution of all the terms can easily be summed and found to give a theory with no phase 
transitions at leading order. The leading term in the free energy of this theory corresponds 
in the string picture to a summation over unbranched covers of a torus by a torus. 

For the spherical case G = 0, however, the situation is more complicated. In the partition 
function, terms corresponding to Young diagrams with n boxes give rise to terms of order 
A^^". The free energy of the string expansion, given by the logarithm of the partition function, 
has a leading term of order A^^. However, an infinite number of string diagrams contribute 
to the coefficient of this leading term, and so a phase transition in the leading order term 
cannot be ruled out. 



The possibility of a phase transition in the genus free energy was investigated in |I4, 



T5| . These papers considered the U{N) theory, using the continuum variables 



to describe Young diagrams in the large N limit. Changing variables to 



h{x) = —n{x) + X — 1/2, 



(2.8) 
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the partition function can be rewritten as a functional integral 

Z{A) = J Vh{x) e-^'^l'^^'M, (2.9) 

where 

S[h{x)] = - [\x dy ln[h{x) - h{y)] + - [\x h{xf - —. (2.10) 
JO JO 2 Jo 24 

This action is familiar from Hermitian matrix model theories, where the log term arises from 
the Van der Monde determinant factor associated with integration over the angular degrees 
of freedom in the matrix. The saddle point of this action is described by the equation 

where uih) = dx/dh. Since x ranges between and 1, u is normalized so that J u = 1. 
Because the variables nj are nonincreasing, u must satisfy the constraint u{h) < 1. 

The integral equation (|2.11| ) has a solution given by the Wigner semicircle law. However, 



for values of A > n'^, the Wigner semicircle solution violates the constraint u < 1. In order 
to deal with this problem, Douglas and Kazakov (DK) adopted the Ansatz that the function 
u{h) should be identically equal to 1 in a region [—b,b], and should be between and 1 in 
the regions {—a, —b) and (6, a) for some < 6 < a. This Ansatz was physically motivated 
by the idea that above the critical area, the density function u would become saturated 
to its maximum value of 1 in a region symmetric around the origin, due to the symmetry 
of the quadratic potential. Using this Ansatz, DK solved the integral equation ( |2.11 ) in 



the nonconstant regions of u, and found an expression for u{h) and the leading term in the 
derivative of the free energy in terms of complete elliptic functions. They found a third order 
phase transition separating the large and small A phases at A = vr^, and verified that the 
free energy in the large area phase agrees with the string expansion. 

In subsequent work ||T^, it was pointed out by Minahan and Polychronakos (MP) that 
the symmetric form of the DK Ansatz effectively restricts the saddle point to the Q = 
U{1) charge sector of the U{N) theory. They considered a more general asymmetric Ansatz, 
in which u has support in the region [d, a) , and is identically 1 in the region [c, b] where 
d < c < b < a. By solving the integral equation ( p.ll|) using this Ansatz, they found a more 
general set of solutions for the U{N) theory, where for a fixed value of A there is a range 
of solutions corresponding to different values of the U{1) charge Q. MP did not ascertain 
whether the value of the free energy was larger or smaller on these additional solutions, so 
that it is not clear whether these other possible saddle points affect the free energy of the 
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theory. However, it was shown in |^ that the free energy of the complete SU (N) theory is 
equal to the free energy of the Q = U (N) theory in the large area phase where the string 
expansion is valid. Thus, the results of DK can be taken to be a complete description of the 
phase structure of the large N SU (N) theory on the sphere. This result seems to indicate 
that in fact, for a fixed A, the DK saddle points represent the extremum of the action with 
respect to Q. 

In this paper we analyze the chiral version of the large N theory using the matrix model 
approach. The precise theory we will investigate is the chiral version of the Q = U{N) 
theory, which corresponds to a string theory of orientation-preserving maps with branch 
points but no tube or handle singularities. To define a matrix model for the chiral theory, 
we choose an Ansatz for which u selects representations with row lengths rij = for all i 
greater than some fixed value k. In the continuum variables, this condition corresponds to 

u{h) = l, foTc<h<l/2, (2.12) 

where u{h) has support on the interval (rf, 1/2) for some d < c < 1/2. The remainder of 
this paper is devoted to using ( |2.12| ) to calculate the saddle point and determine the phase 
structure of the chiral theory. 



3 Single-cut Solution 



The simplest Ansatz for a function u{h) satisfying (|2.12| ) is that m < 1 in the region {d, c) 



This will result in a solution that is the analog of the Wigner semicircle saddle point solution 
in the low area phase of the coupled theory 0. Surprisingly we will find that it describes 
both the small and large area phases of the chiral theory. This Ansatz gives rise to the 
integral equation 

^+,„izV! = prd.i^(fl. (3.1) 

2 h — c Jd h — s 



The general approach to solving an integral equation of this type [|T^ is to define a function 



f{h)=j\s^. (3.2) 
Jd h — s 

Note that / is analytic on the complex h plane, with a cut along the interval [d, c] . Defining 
f+{h) and f~{h) to be the limiting values of / as one approaches the cut from above and 
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below respectively, we have 



and 



(3.3) 



(3.4) 



Let g{h) be a function which is valued on the cut plane, and let g+{h) and g-{h) be the 
limiting values of g{h) as one approaches the cut from above and from below respectively. 
We choose a g which satisfies g^ + g^ = 0. It follows that 

u{h) 



m- 



[if/gU + if/gU 



and 



But this indicates that 



g+{h) 2 

if/gU-if/g). /+ + /_ 



27ri 

g(h) 



ds 



2mg+ 

ni {h - s)g+{s)' 



-1^ 

m g+' 



For our problem, a simple choice for a function g satisfying + gi_ = is 



and thus 



fih) 



f ds 



gM. 

2m I {h — s)g{s) 



g{h) = J{h-c){h-d), 



±:^{h-c){h-d)fds 



s-l/2 
s—c 



{h — s)J{s — c) (s — d) 



(3.5) 

(3.6) 
(3.7) 

(3.8) 
(3.9) 



where the integral is taken around the cut (c, d). Deforming the contour, we get contributions 
from the poles at h and oo and the integral around the log cuts. 



Ah 



2 ■ h-c 
We can now calculate u{h) in the interval {d, c) 



1/2 J{h -c){h- d) 

ds . 

" {h - s)J{s - c){s - d) 



(3.10) 



u{h) 



1 

71 

1 

TT 



A^{c-h){h-d) 



r.l/2 



'c — h){h — d) 



+ tan" 



{h — s)^^{s — c)(s — d) 

^ ( ^{c-h){h-d){l-2c){l-2dy 
h- {h + l/2){c + d) + 2cd 



(3.11) 
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Figure 1: single-cut solutions 



where we take values of tan~^ to be between and vr. 

The variables c and d are of course not arbitrary. Since u is normalized to J u{h) dh = 1, 
we know that for large h, f{h) satisfies 



1 - 1/2 + c 
h 



+ (■ ■ ■)h-' + {-2F'{A) + 1/24 + cy3)h-'^ + 0{h 



(3.12) 



where we note that F{A) is the free energy, F{A) = ln{Z {A)) / N"^ . Expanding ( |3.10[ ) in 
powers of h and performing the necessary integrations, we find that c and d satisfy the 
conditions 

A 



[c + d) = In 



[l-2c){l-2d) + l-d-c 



1 



A^^^ + ^{l/2-c){l/2-d). 



(3.13) 



(3.14) 



In Figure |^ we have graphed the solutions to (|3.13|) and (|3.14|) on the c-A plane. These 
solutions lie on two disconnected curves. We observe that in the large area regime, as A —>■ oo, 
we have c,d ^ ~V2- As A 0, on the other hand, we have c — 1/2, d — oo. In section 
4 we show that the large A solution to these equations reproduces the string solution. 
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For a solution of the constraint equations to give a physically acceptable saddle point we 
must further verify that u{h) < 1. Unphysical solutions are denoted in Figure |l| by dashed 
lines. Checking the solutions in the large A region, we find that for sufficiently large c, the 
function u is larger than 1 for h c—. The A value for which this first occurs is the point 
where dw/d/i = at /i = c. Furthermore, it is straightforward to verify that this is precisely 
the point where dA/dc = drf/dc = 0. This implies that c and d satisfy the relation 



1/2 -c 



+ 3c + rf-2 = 0. 



(3.15) 



^ 1/2-ci 

Combining this with ( p.l3| ) and ( |3.14|) we find that the area at this point is approximately 



A+ ^ 10.189 



(3.16) 



Similarly, we find that for sufficiently large A, the solutions on the small A curve in 
Figure are unphysical because the function u develops a "bump", where u{h) > 1. The 
boundary of the physical region occurs at the point where the equations 



u{e) = 1 u'{e) = 



(3.17) 



have a simultaneous solution for some d < e < c. These equations simplify to the simulta- 
neous equations 



cos 



'A^J{c-e){e - dy 

^(l/2-c)(l/2-rf) 
1/2 - e 



(e-d)(l/2-c) 
\\{c-d){l/2-e) 



^{2e-c-d) 



(3.18) 



(3.19) 



Again combining with ( p. 13 ) and ( 3.14 ), the solution to these equations corresponds to an 
area of approximately 

A_ ^ 9.426 (3.20) 

In Figure H, we have graphed the function u{h) at the two values of area Aj^ and A_. These 
represent the boundaries of the physical region of the single-cut solution. 

Given the solutions we have found to the integral equation (|3.1|) which satisfy the con- 
straint u{h) < 1, we can use ( |3.10| ) and ( |3.12D to calculate F'{A), the derivative of the free 




Figure 2: Function u{h) at critical points 
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energy with respect to the area, for the solutions of the single-cut Ansatz. We have 

{c-df{5c'^ + 6cd + 5d^)A 



- 2F'{A) 



-1/12 



+ 



1/2 ^3 

ds— 



256 

s\c + d)/2 - s{c - df/S -{c + d){c- df/lQ 



(3.21) 



(s — c)(s — d) 



Performing the integration and using ( 3.14|) , we have 



F'iA) 



A{c-dy^ 



{1 + c + d 



9c2 + Ucd + 



,c + d ^ 3c2 + 2cc/ + 3rf2^ 



(3.22) 



384 ' 4 ' ' 24 48 ' 

We have thus found that for any area outside the region {A_,A^) there is a unique 



solution to the integral equation (|3.1| ) which satisfies u{h) < 1 for d < h < c. We have 
calculated the derivative of the free energy for these solutions. However, for A_ < A < Aj^ 
we have found no solution corresponding to a function u representing a master field for the 
large Young diagram. In Section 5, we will show that allowing u to be identically 1 in 
another interval provides a solution in the intermediate area regime. 



4 Enumeration of Sphere Maps 

Now that we have an analytic expression for the derivative of the free energy of the chiral 
theory for A > Aj^ and for A < A_, we expect that in the large area regime this formula 
should agree with the free energy of the chiral QCD2 string. The chiral string theory has a 
free energy which is expressed as a double expansion in e""^/^ and A. Specifically, the free 



energy of the chiral theory can be written ||Tl 



Y.MA)x\ (4.1) 



n=l 

where 

x = e"?, (4.2) 
and the functions 0„(A) are polynomials of degree 2n — 2 in A. To check the correspondence 



A 



with (|3.22|) we must expand our solutions for c, d, and F' as double expansions around 
A = 00. Writing 

c=-l/2 + e d = -l/2-5, (4.3) 
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the conditions ( p. 13 , 3.14 ) on c, d become 



^-A/A^{e-5)A/A ^ 

2(1 + J(l-e)(l + 5)) + 5-e 



and 

\2 



l-V(l-e)(l+i) = ^(i±^L. (4.5) 

We can rewrite these two equations as 



e + 5 = V^(A- ^il±3l + s-e^ e(-^)^/^ (4.6) 

These equations can be solved order by order in ^/x, giving c and d as power series' in x^/^ 
with the coefficients of x"/^ being polynomials in A of order n — 1. Computing up to n = 2, 
for example, we have 

c = -l/2 + 2x^^^ + {A-2)x + {A'^-AA + 2)x''/^ 

+ {A' - QA^ + 8A- 2)x^ + 0{x^'^) 

d = -l/2-2x^/^ + {A-2)x - {A^ -4:A + 2)x^/^ (4.8) 
+ {A^ - QA^ + 8A- 2)x^ + 0{x^'^). 



Plugging these expansions into ( p.22| ), we have an expansion for F'[A) in x and A. The 



lowest order terms in this expansion give 

F\A) = -\x-{^-^-^ + \)x' + ... (4.9) 

We have calculated this expansion out to order x^^, and find exact agreement with the 
coefficients in the polynomials 0„ of the string expansion. The polynomials 0„(A) for n < 7 
are given explicitly in 0]. 

Thus, we have verified that in the large area regime, the matrix model we have constructed 
correctly reproduces the string theory of chiral QCD2. We now have an exact formula for 
the free energy of this theory, which we can use to calculate to arbitrary order the terms 
in the expansion. A particularly interesting calculation is the determination of the leading 
coefficients in the polynomials (f)n{A). For a fixed value of n, the polynomial (pn is of degree 
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2n — 2, and the coefficient 7„ of y4^"~^/(2n — 2)! is equal to the number of topologically 
distinct holomorphic maps from 5*^ to 5*^ with 2n — 2 elementary branch point singularities 
whose images are fixed. It was conjectured by Gross and Taylor [11| that these coefficients 
are given by 

n"-3(2n-2)! , 

7n = 4.10 

n! 

We can now use the exact expression for the free energy (|3.22|) to prove that these are 
precisely the leading coefficients in the matrix model free energy expansion, and thus by 
implication that the coefficients 7„ of the QCD2 string are indeed given by (|4.1CI|) . 

To calculate the desired coefficients in the free energy expansion, it will be useful to 
consider a simplification of the theory described above in which subleading terms in the 
polynomial coefficients of powers of x are dropped. We denote the free energy of the simplified 
theory by F. This free energy has a derivative which is of the form (continuing to drop 
subleading terms) 

n 

In the string theory language, this simplification corresponds to removing ll-points from the 
theory. From the form of equations (|4.6|) and ([4.7|) , it is straightforward to show that the 
expansions of e and 6 are of the form 

00 

e = ^x"/2(e„A"-i + C(A"-2)) (4.12) 

n=l 

00 

S = ^a;"/2(c;^A'^~i + C)(A"-2)). (4.13) 

n=l 

Rewriting ( |3.22|) in terms of e and 6, and dropping terms which are of lower order in A for 
each power of x, we have 

^ ^ 2(i-sr-(i+sr _ 1 1 

V ; 8^ ' 2A^ ' ^ ^ 



where 



e„A"-V"/2 (4.15) 

n=l 

5] rf„v4"- V/l (4.16) 



n=l 



Dropping lower order terms in equations ( |4.6| ) and ( [4.7| ), we have 

i + S = 4x'/^e^'-')^/' (4.17) 
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and 

^-i^^. (4.18) 

Combining these equations, 

^ = ^^l_3. = A2^e«. (4.19) 
We will now find it useful to recall some properties of the generating function 

A famous theorem due to Cayley [|l^ asserts that c„ = n"~^ is the number of distinct 



connected graphs with n labeled vertices and n — 1 edges. If we remove a fixed vertex from 
every graph of this type, and sum over all possible graphs by summing over the number k of 
edges connected to that vertex and over all possible subgraphs connected to each edge, we 
can derive the recursive equation 

Cn = E E fc'iT'T.l,, ]!^-^-^- (4-21) 

This equation implies that f{z) satisfies 

f(z) = ze^^'^ (4.22) 

If we were to label edges instead of vertices, the number of distinct graphs would be c„/n 
(where graphs with symmetry are counted with a weight of the inverse of the order of 
symmetry). By removing a fixed edge from each graph, and counting all graphs connected 
to that edge as above, we arrive at the equation 

"■"^ (n-2V 
cJn^l/2Y. ^_^,^J,_,^, Cc.-, (4.23) 

This gives rise to the equation for the generating function, 

oo 



/W' = 2E7 (4-24) 

^1 (^-2)! 



From ( ^^221) and (|4l9| ), we see that 

^ = f{A'x). (4.25) 
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We can now rewrite the derivative of the free energy in the form 



P,A).i^^-^.-E^A^"-^'. ,4.26) 

This is precisely the desired result, and proves that the coefficients in ( [4.10|) describe the free 
energy for the large area phase of the matrix model described in the previous section. The 
equivalence of this matrix model with the chiral string expansion of QCD2 indicates that the 
number of topologically distinct degree n covers of the sphere by a sphere with 2n — 2 fixed 
branch points is exactly n"~^(2n — 2)\/n\. This result was used (but not proven) in |[ll[] to 
argue that in the strong coupling string expansion the entropy of branch point singularities 
drives the phase transition on the sphere, both in the chiral and coupled theories. 



n-3 



5 The Double-Cut Solution 

In Section 3, we derived a single-cut solution to the equation of motion for both small 
{A < A^) and large areas {A > A^). For intermediate areas A_ < A < A^ there are 
no physically viable single-cut solutions. This is due to the fact that the functions u{h) 
describing the single-cut solution for the areas in this range are "unphysical" ; they do not 
satisfy the constraint u{h) < 1 and thus cannot correspond to Young diagrams. 



This situation is analogous to that which arises at the DK phase transition point j^, |15 
In the coupled theory, the simple saddle point solution described by the Wigner semicircle 
law gives a good description of the theory for areas less than the critical area A = n^. 
However, above the critical area the Wigner semicircle violates the constraint u < 1, and 
does not correspond to an acceptable solution. In the language of p|, |T^, this theory is 
described as a system of fermions moving in 1 dimension in an external quadratic potential 
and interacting with a logarithmic potential. When the critical area is reached, the fermion 
density is saturated at /i = and condensation occurs giving a function u{h) which is 
identically 1 in a small region. 

In terms of the fermion language, we are here considering a system where the potential 
is again quadratic but where there is an impenetrable barrier at the point h = 1/2, at which 
fermions tend to congregate. As the area is increased, the strength of the external quadratic 
potential increases and the particles move towards the origin h = 0, creating a "bump" in 
u{h). Eventually, the first critical point A- is reached, and the fermion density reaches unity 
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at a point h = e. Just as the full theory developed a condensate where u = 1 beyond the 
critical point, we might expect a similar effect in the chiral theory above A = A_. However, 
unlike the situation which Douglas and Kazakov considered, our system is not symmetric, so 
we cannot expect any particular relation between the endpoints of the condensation region. 
Thus, we will adopt an Ansatz for the solution of ( p.ll| ) in which u{h) has support on the 
region (rf, 1/2), with 



u{h) = 1 for c<h<b or a<h<l/2 



(5.1) 



for some values d < c < b < a. Note that u{h) is subject to the constraint u{h) < 1 
throughout the regions {d,c) and (6, a). 

A solution, u{h), satisfying the double-cut Ansatz ( [5.1D may be found by adapting the 
method described in Section 2. The integral equation for the saddle point thus reads. 



Ah {h-l/2){h-b) 
{h~a){h-c) 



P / ds 



h 



(5.2) 



where the region of integration is R = {d,c) U (6, a). As before, it is useful to define the 
analytic function 



fih) 



ds 



MS 



h 



(5.3) 



and, letting g{h) = ^{h — a){h — b){h — c){h — d) as appropriate for the double-cut Ansatz 
( p.l| ), we find that we may rewrite f{h) as 



1/2 



ds 



{h - s)g{s) 



+ 9{h) 



ds 



{h - s)g{s) 



(5.4) 



where we have deformed the contours giving terms corresponding to the pole at s = h and 
to the logarithmic discontinuities on the intervals (a, 1/2) and (c, 6). We may express f{h) 
from ( ^.4| ) in terms of elliptic functions. This is particularly useful for writing u{h) in the 
region R. We find 



u{h) = — 



2 

np 



\ {a-h){h-b) 



+ 



{a-h){h-b) 
\\ {h-c){h-d) 



{b-d){h-a) 



(5,5) 



where 



sm 



{b- d){l/2- a). 
V (a-c/)(l/2-6)- 



{b — c)(a — d) 
^ {a-c){b-d) 



p= J(^a-c){b-d) (5.6) 
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and where the notation for the elhptic functions is that used in Byrd and Freidman WE 



Thus F,E and 11 are the elhptic integrals of the first, second and third kind, respectively, 
and when the angle u does not appear as one of the arguments of an elliptic function the 
symbol then represents the associated complete elliptic function (note K{q) is the complete 
elliptic function associated to F). 

As described in Section 3, the parameters a, b, c, d are not arbitrary but are subject to 



conditions related to the definition of f{h) in (|5.3| ). As was done in Section 3, we study the 
asymptotics of / to determine these conditions on a,b,c,d. Expanding f{h) about h ^ oo 
and isolating the terms of order 0{h), 0{1) and 0{h~^) respectively we find, 

A = -iKiq)-F{u,q)) (5.7) 
P 

{a - b)U{u, q) + {b- W/2)F{u, q) = [c - dM^, q) + {d - W/2)K{q) (5.8) 



'l/2-a)(l/2-c)(l/2-d) , (a-6-c + rf)^ 

' ' ' +p{E{q)-E{u,q))+^- L{^K{q) - F{v,q)) 



\ 1/2 -b ' '^^ ^ ' 4p 

(5.9) 

where W = a + b + c + d. Note that since there are four unknowns a, 6, c, d, for fixed area A 
there will be at most a one parameter family of solutions to these three equations. Also note 
that as one takes a b [or c ^ b) two of the three equations above become the constraint 
equations (|3.13|) , (|3.14|) of the single-cut solution. 

The solutions to this double-cut problem which satisfy the constraint u{h) < 1 for all 
h are depicted in Figure |^. We find that there is a continuum of double-cut solutions (the 
shaded region) for each A > A_. The upper boundary of the region of solutions is given by 
the curve describing the single-cut solution above A+. The cusp point at the bottom of the 
region is precisely the lower critical point. Note that there the value of c is the value of h 
at which the "bump" in u (see Figure ^ touches u{c) = 1, and not the value defined in the 
single-cut solution, which is now taken by a at that point. 

The constraint u{h) < 1 must be satisfied for all h in order for the solution u{h) to cor- 
respond to a Young diagram at large N. This condition can be used to describe analytically 
the remaining boundary curves of the physically acceptable region in Figure |^. As one varies 
a,b,c,d in such a way as to satisfy (|5^ ), (|5.8|) , and ( ^.9|) , the solution (|5.5| ) first fails this 
condition at the points b and c. Thus, we can find the boundaries of the double-cut problem 
simply by studying the slope of the solution u{h) at the points b and c. The boundary of 
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-0.4 
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-0.1 



0.1 



Figure 3: double-cut solutions 



the double-cut solutions arising from the condition that u{b-\-) < 1 is 
(a - b)iKiq) - F{u, q)) = {a - c){E{q) - q)). 



(5.10) 



This is the boundary curve at the bottom of the physical region. An analysis of the slope 
of the solution u{h) near the point c, {i.e., requiring u{c—) < 1) yields an equation for the 
final boundary of the double-cut region. It reads 



p{b 



'l/2-d){l/2-a) 
\ (l/2-c)(l/2-6) 



(a _ c)(c - d){K{q) - F(z/, q)) - p\E{q) - E{u, q)). {5.11] 



One finds that the solutions at these boundaries are indeed legal solutions, so the double-cut 
region is closed. It is easy to verify that the u{c—) < 1 boundary runs precisely between the 
y4+ and A_. 

Having found these new saddle points of the action ( |2.10| ) we now study their contribution 
to the partition function. To begin with, we can see from Figure |^ that for all A > there 
is indeed a one parameter family of solutions to the saddle point equation. However, not all 
these solutions are actually saddle points, since the saddle point equation is only satisfied in 
the interior of the regions {d, c) and {b, a). To determine which of the values is the physical 
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saddle point for any given A it would be necessary to compute the exact free energy F and 
to maximize this quantity across all solutions corresponding to that fixed value of A. 

Unfortunately, however, an exact calculation of the free energy is technically difficult. 
The situation here is analogous to that which arises in the full QCD2 theory. In that case, it 
was shown by Minahan and Polychronakos |15| that for every value of A above the critical 



value 71^ there is a family of solutions to the saddle point equation, parameterized by the 
parameter Q corresponding to the U{1) charge, which is essentially the first moment of u{h). 
In that case also, a direct calculation of the free energy is difficult. Symmetry arguments 
would indicate that the DK solutions with U{1) charge Q = are the physical saddle points 
in this theory. The correspondence of the strong coupling expansion of the DK solutions 
with the QCD2 string gives further weight to this conclusion; however this result has not 
been rigorously proven. 

In the chiral theory, we no longer have a reflection symmetry to rely on; however, the 
correspondence of the strong coupling expansion of the single-cut solution with the chiral 
string theory is a strong indication that the free energy is extremized by the single-cut 
solutions above A = A^. We have numerically integrated the free energy in that region 
and found that indeed F seems to be maximized on the single-cut solution for any flxed 
A > y4+. The behavior of the theory in the intermediate phase, however, presents a much 
more subtle problem. Numerically integrating F indicates that the physical solutions in the 
intermediate phase lie near or on the boundary ( p. 111 ). The numerical analysis of F in this 



region is extremely delicate, presumably due to the existence of phase transitions at A_, A^. 

Thus, although we have fairly solid evidence that the single-cut solution describes the 
physical saddle point above A+, we are unable to give an exact formula for the curve corre- 
sponding to the physical solutions connecting the two critical points. Nonetheless, we can 
compute the derivative of the free energy along the physical curve using the same method 
as was used in Section 3 to calculate F'{A) in the single-cut region. From the fact that the 
physical saddle point obeys the classical equations of motion, we flnd that the derivative of 
the free energy is related to the second moment of the density function u{h) by 

F' (A) = J dhu{h)h^ + 1/24. (5.12) 

After some algebra one ascertains that 

/ dhu{h)h^ 
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[1/2 -m + w)^^^^ 



6 ^\ 1/2-6 



(l/2-a)(l/2-c)(l/2-c/) _ , 



+ (3W^2^2 ^ ei^XFZ + ^ 2X2^2 + 2X^Z^ + x4)-^^M_ZiM) (5.13) 

where 

W = a + b + c + d X = a-b- c + d Y = a + b- c- d Z = a-b + c- d (5.14) 
and where 

M = {3W^ + X^ + + Z'^)/A8 (5.15) 



Using the constraint equations ( p.7| , |5.8| , |5.9| ) allows us to write this expression in terms of 



algebraic functions of the area A and the points a, b, c, d. 

{l/2-b){l + W) 



dhu{h)h' 



(l/2-a)(l/2-c)(l/2-rf) , 
6 \ 1/2-6 

+ {6WXZY + Y^Z^ + X^Y^ + X^Z^)A/7Q8 (5.16) 



As a simple check note that in the limit b a this expression leads to the result of (|3.22| ) 



We have thus given an explicit formula for the derivative of the free energy in the double-cut 
region which is applicable to whatever curve describes the physical solutions in that region, 
although we have been unable to describe the curve analytically. If one could derive an 
analytic formula for the physical solutions in the intermediate region, ( ^.16| ) could be used 
to calculate the order of the phase transition at the two critical points. 



6 Conclusions 

We have developed a matrix model that describes a single chiral sector of the QCD2 string. 
The analysis indicates that there are three distinct regimes for solutions of the saddle point 
equation of the chiral theory on the sphere. We have shown that there exists a simple 
single-cut solution which exists for large and small areas and which agrees with the string 
picture in an expansion around large area. However, this single-cut solution ends (stops 
leading to physically viable solutions) as one decreases the area from infinity and also as 
one increases the area from zero. Thus the single-cut solution cannot describe the system 
for areas A_ < A < We have shown that for these intermediate saddle point 



21 



may be found as a double-cut solution to the (integral) equation of motion. This double-cut 
solution limits to the single-cut solution in the large area regime and connects the upper and 
lower branches of the single-cut solution. For any fixed A > there is a one-parameter 
continuum of saddle point solutions. When A > A+, we have fairly conclusive evidence 
that the single-cut solution dominates the theory. However, in the intermediate regime we 
cannot calculate exactly where the curve describing the physical saddle point lies. In fact, 
although we believe there is only a single curve with maximum free energy connecting the 
two points A^ and A_, it is conceivable that there may be multiple physical saddle points 
in this intermediate region, although numerical evidence indicates that there is indeed only 
a single such curve, which lies near or on the u{c—) < 1 boundary ( ^.11|) . 

Despite our uncertainty regarding the exact location of the intermediate saddle point 
curve, our analysis gives a fairly complete picture of the theory at any A. The dominant 
stationary solutions are the single-cut solutions for A > A+ and A < A^ and there are likely 
to be a simple trajectory of solutions in the double-cut region that interpolate between these 
single-cut solutions. The fact that the physical curve of solutions in (c. A) space, as well 
as the solutions u{h), undergo discontinuities at the points A^ and A_ indicates that there 
are probably phase transitions at these points. The behavior of the density function u{h) 
at the point A^ is closely analogous to that of the full QCD2 density function at the DK 
phase transition point, suggesting that the chiral theory probably has a third order phase 
transition at this point. However, an exact characterization of phase transitions in the chiral 
theory will not be possible until a method is found for determining which of the solutions in 
the double-cut region actually dominate the partition sum. If this could be accomplished, it 
would then be straightforward in principle to compute the order of the phase transition at 
the critical points simply by comparing power series expansions of the free energy F{A) on 
either side of the critical points. 

The single-cut solution of the matrix model gave us an exact analytic formula for the 
large area regime of the chiral theory. Because this exact formula has a series expansion 
which is equal to the leading (A^^) term in the free energy of the chiral string theory, we were 
able to relate the terms in this expansion to sums over maps between Riemann surfaces. 
In particular, a calculation of the leading terms in the area polynomials of the expansion 
allowed us to derive a previously conjectured result on the number of sphere to sphere maps. 
An interesting direction for further work might be to find further structure in the analytic 
equations for the large area chiral theory and to relate this structure to the geometry of 
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string maps. 

It is remarkable that the location of the critical point bounding the large area phase is 
relatively insensitive to the details of the theory in question. The result proven here for 



the number of sphere to sphere maps was used in [|TT| to show that in a chiral theory with 
branch points, but without i7-points, the critical point occurs near A = 11.9. In this paper 
we found that the i7-points in the chiral theory shift the critical point to near A = 10.2, 
which is extremely close to the DK critical point in the coupled theory at A = n'^ ^ 9.9. 

One particularly surprising feature of the results given here is the fact that although the 
chiral theory apparently has three distinct phases, the large and small area phases are in fact 
described by the same analytic expressions, indicating some type of duality between these 
two phases. This situation also bears an striking similarity to the results of [Tl| where the 
convergence properties of the string expansion were studied for both the chiral and coupled 
QCD2 theories. In that analysis, it was shown that the string expansion of the chiral theory 
has a radius of convergence approximately equal to the point A^, but that the string series 
also converges for very small areas. A numerical comparison, however, shows that the results 
of the matrix model disagree with the explicit string sum for small area. The small area 
phase we have described in this work of course lies outside the radius of convergence of 
the large area phase. Thus, the fact that the string sum did converge for small A seems 
mysterious. 

In order to further understand the chiral QCD2 matrix model we have described here, it 
would be interesting to study directly the free energy of this model. A way of analytically 
describing the dominant saddle points for A_ < A < A^ in this model would be invaluable 
for the computation of the order of the phase transition across the critical loci. Furthermore, 
an understanding of how to compute F directly might shed some light on the significance of 



the one-parameter families of solutions described in |15| for the full QCD2 large N theory. 
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